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Summary 

A set of compressible flow relations for a ther- 
mally perfect, calorically imperfect gas is derived for 
a value of c p , expressed as a polynomial function of 
temperature, and developed into the thermally per- 
fect gas (TPG) computer code. The code produces 
tables of compressible flow properties similar to those 
in NACA Rep. 1135. Unlike the tables of NACA 
Rep. 1135 that are valid only in the calorically per- 
fect temperature regime, the TPG code results are 
also valid in the thermally perfect, calorically imper- 
fect temperature regime, which considerably extends 
the range of temperature application. Accuracy of 
the TPG code in the calorically perfect temperature 
regime is verified by comparisons with the tables of 
NACA Rep. 1135. The small differences observed in 
this comparison result from inaccuracies of the tab- 
ulated values of NACA Rep. 1135, which originate 
from the calorically perfect gas assumption. In the 
thermally perfect, calorically imperfect temperature 
regime, the TPG code is validated by comparisons 
with results obtained with the method of NACA 
Rep. 1135 for calculating the thermally perfect, calor- 
ically imperfect compressible flow properties. The 
advantage of the TPG code is its applicability to any 
type of gas (monatomic, diatomic, triatoinic, poly- 
atomic, or any specified mixture thereof), whereas 
the method of NACA Rep. 1135 is restricted to only 
diatomic gases. 

The temperature limits for application of the 
TPG code are also examined. The thermally perfect 
behavior of any gas species will terminate at some 
combination of low temperature and high pressure 
because of the presence of real gas effects and at some 
value of high temperature because of dissociat ion (ex- 
cept monatomic gases). Beyond these temperature 
limits, the TPG code may no longer accurately pre- 
dict gas properties. However, within the temperature 
range in which the thermally perfect assumption is 
valid, the TPG code can accurately predict the com- 
pressible flow gas properties of any given gas species 
or mixture of species. 

Introduction 

The traditional computation of the one- 
dimensional (1-D) compressible flow gas properties 
is performed with one-dimensional calorically per- 
fect gas equations such as those of NACA Rep. 1135. 
(See ref. 1 .) If the gas of interest is air, then the ta- 
bles of compressible flow values provided in NACA 
Rep. 1135 can be used. These tables were gener- 
ated from the calorically perfect gas equations with 
a value of 1.40 for the ratio of specific heats 7 . The 
application of these equations and tables is limited 


to that range of temperature for which the calorically 
perfect gas assumption is valid. However, many aero- 
nautical engineering calculations extend beyond the 
temperature limits of the calorically perfect gas as- 
sumption, and the application of the tables of NACA 
Rep. 1135 can result in significant errors. The tem- 
perature range limitation is greatly minimized by the 
assumption of a thermally perfect, calorically imper- 
fect gas in the development of the compressible flow 
relations. (For simplicity, the term thermally per- 
fect will be used throughout this paper to denote a 
thermally perfect, calorically imperfect gas.) This 
paper presents a computer code which implements 
one-dimensional compressible flow relations derived 
for a thermally perfect gas. 

A calorically perfect gas is by definition a gas for 
which the values of specific heat at constant pres- 
sure c p and specific heat at constant volume c v are 
constants. Therefore, in the derivation of the com- 
pressible flow relations for a calorically perfect gas, 
the value of c p was assumed constant. The calorically 
perfect gas equations are derived in many compress- 
ible flow textbooks and are summarized in NACA 
Rep. 1135. The accuracy of these equations is only as 
good as the assumption of a constant c p (and there- 
fore a constant 7 ). For any diatomic or polyatomic 
gas, the value of c p actually varies with temperature 
and can only be approximated as a constant for a rel- 
atively narrow temperature range. At some point as 
the temperature increases, the c p value begins to in- 
crease appreciably due to the excitation of the vibra- 
tional energy of the molecules. This phenomena oc- 
curs around 450 to 500 K for air. Thus, above 500 K, 
the use of 7 = 1.40 for air in the calorically perfect 
gas equations will yield noticeably incorrect results. 
The variation of c p with temperature (and only with 
temperature) continues up to approximately 1500 K. 
In this temperature range, where the value of c p is 
only a function of temperature, air is considered a 
thermally perfect gas. Therefore, the derivation of 
the compressible flow relations for a thermally per- 
fect gas must account for a variable heat capacity. 
Above 1500 K for air and at some relatively high 
temperature for all diatomic and polyatomic gases, 
dissociation of the molecules starts to occur. When 
dissociation commences, the value of c p becomes a 
function of both temperature and pressure and the 
gas is no longer considered thermally perfect. The 
calorically perfect and thermally perfect temperature 
limits cited for air illustrate the much wider temper- 
ature range of application available from thermally 
perfect compressible flow relations. 

The NACA Rep. 1135 presents one method for 
computation of the one-dimensional compressible 



flow properties of a thermally perfect gas. In this 
approach, the variation of heat capacity due to the 
contribution from the vibrational energy mode of the 
molecule is determined from quantum mechanical 
considerations by the assumption of a simple har- 
monic vibrator model of a diatomic molecule. With 
this assumption, the vibrational contribution to the 
heat capacity of a diatomic gas takes the form of 
(ref. 2) 


(1L - ( 

where 0 = hvjk , h is Planck’s constant, v is the 
characteristic frequency of molecular vibration, and 
k is the Boltzmann constant. The complete set of 
thermally perfect compressible flow relations based 
upon the Cp variation given in equation (1) is listed 
in NACA Rep. 1135 “Imperfect-Gas Effects.” Tables 
of these thermally perfect gas properties are not 
provided because each value of total temperature 
Tt would yield a unique table of gas properties. 
The NACA Rep. 1135 does provide charts of the 
thermally perfect air properties normalized by the 
calorically perfect air values and plotted versus Mach 
number for select values of total temperature. This 
approach can only be applied to diatomic gases (e.g., 
N2, O2, and H2) because of the functional form 
which describes the variation of heat capacity with 
temperature. 
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temperature and is presented in “Derivation of Ther- 
mally Perfect Gas Equations” in this report. This set 
of equations was then coded into a computer program 
referred to as the thermally perfect gas (TPG) code, 
which represents the end product of this research ef- 
fort. The output tables of the TPG program are 
structured to resemble the tables of one-dimensional 
compressible flow properties that appear in NACA 
Rep. 1135; the difference between the two sets of ta- 
bles is that the TPG code is valid in the thermally 
perfect temperature regime as well as in the calor- 
ically perfect temperature regime. This code serves 
the function of the tables of NACA Rep. 1135 for any 
gas species or mixture of gas species (e.g., air) and 
significantly increases the range of valid temperature 
application. 

Symbols 

a 
A 


c v 

c v 

e 

h 


speed of sound 

cross-sectional area of stream tube 
or channel 

coefficients of polynomial curve fit 

specific heat at constant pressure 
specific heat at constant volume 
internal energy per unit mass 

enthalpy per unit mass, e -f pv; 
Planck’s constant 


Because equation (1) is applicable to only di- 
atomic gases, a different method for computing the 
one-dimensional compressible flow values of a ther- 
mally perfect gas was developed and is described in 
this paper. This method utilizes a polynomial curve 
fit of Cp versus temperature to describe the variation 
of heat capacity for a gas. The data required to gen- 
erate this curve fit for a given gas can be found in tab- 
ulated form in several published sources such as the 
NBS “Tables of Thermal Properties of Gases” and 
the “JANAF Thermochemical Tables.” (See refs. 3 
and 4, respectively.) Actual coefficients for specific 
types of polynomial curve fits are also published in 
NASP TM-1107 and NASA SP-3001. (See refs. 5 
and 6, respectively.) These curve fits, which are 
based upon tables of standard properties of gases, 
enable the application of this method to any type 
of gas; monotomic, diatomic, and polyatomic (e.g., 
H2O, C0 2} and CF4) gases, or mixtures thereof are 
included, provided that known data sets of c p versus 
temperature exist for the gases of interest. 

A set of thermally perfect gas equations is de- 
rived for the specific heat as a polynomial function of 


k Boltzmann constant 

V 

M Mach number, — 

a 

p static pressure 

, . pV 2 

q dynamic pressure, ; 

heat added per unit mass 

R specific gas constant 

1Z universal gas constant 

s entropy per unit mass 

T static temperature 

u velocity component parallel to 

free-stream flow direction 

v specific volume, - 

P 

V speed of flow 

w external work performed per unit 

mass 

W molecular weight 


2 


X mole fraction 

Y mass fraction 

Z compressibility factor 

0 ^/\M^X\ 

S increment indicator 

7 ratio of specific heats, — 

© constant of molecular vibrational 

energy 

v characteristic frequency of molecu- 

lar vibration 

p static mass density 

Subscripts: 

1 conditions just upstream of shock 
wave 

2 conditions just downstream of shock 

wave 

t total (stagnation) conditions 

* critical (sonic) conditions 

i ith component gas species of 

mixture 

j j th coefficient of polynomial curve 

fit for Cp 


jmax 

maximum j value 

lim 

temperature limit (upper or lower) 
of polynomial curve fit for Cp 

mix 

gas mixture 

n 

total number of gas species that 
comprise gas mixture 

rev 

reversible 

perf 

quantity evaluated for gas that 
is both thermally and calorically 
perfect 

therm perf 

quantity evaluated for gas that is 
thermally perfect but calorically 
imperfect 

vib 

vibrational contribution 

Abbreviations: 

l-D 

one-dimensional 

JANAF 

Joint Army-Navy- Air Force 

NACA 

National Advisory Committee for 
Aeronautics 

NASP 

National Aero-Space Plane 

NBS 

National Bureau of Standards 

TPG 

thermally perfect gas 


Derivation of Thermally Perfect Gas Equations 
Polynomial Curve Fit for c p 

The selection of a suitable curve fit function for Cp is the starting point for the development of thermally 
perfect compressible flow relations. The form chosen for the TPG code was the eight-term, fifth-order 
polynomial expression given below in which the value of Cp has been nondimensionalized by the species gas 
constant. 

\ = M (^ 2 ) + 42 (^) + 4? + MT)+ a 5 (t 2 ) + a 6 (t 3 ) + a 7 (t*) + a 8 (t 5 ) = Y l AjTi- 3 ( 2 ) 

3 = 1 

For a given temperature range, A\ to Ag are the coefficients of the curve fit. This is the functional form for Cp 
of reference 5, which provides values of A\ to Ag for 15 different gas species. Reference 6 provides similar c p /R 
coefficient data for over 200 different gas species for a five-term, fourth-order polynomial curve fit equation (^3 
to A 7 ). Equation ( 2 ) is also valid with these coefficients provided that A\, A 2 , and Ag are set equal to zero. 
If the pertinent Cp / R versus temperature data fits some other algebraic expression, then that expression could 
easily be substituted for equation ( 2 ) in the TPG code. This issue will be addressed further in “TPG Code 
Description.” The only requirement is that closed-form solutions to both f CpdT and / dT must be known; 
the reason will become apparent as the gas expansion relations are derived for a thermally perfect gas. 
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The TPG code can compute the thermally perfect gas properties for not only individual gas species but also 
for mixtures of individual gas species (e.g., air). This code capability is achieved by calculating the variation 
of the heat capacity for the specified gas mixture. The Cp of a mixture of gases is computed from 

n 

c p,mix = ^ ^ Yi c p,i (^) 

i= 1 

where Y{ is the mass fraction of the ith gas species. The value of Cpj is determined from equation (2) for each 
component species. The following example of standard air, consisting of the four major component gas species 
N 2 , O 2 , Ar, and CO 2 , illustrates how this equation is actually implemented in the code. 

Cp,air = ^N2^P?N2 ^Ar^p,Ar "h ^C 02 ^p,CO 2 (^) 

Equation (4) can be fully expanded with the use of equation (2) to give 

Cj^air ~ ^N2^N2 y ^l,N2 ^ 2 ^ ^N2^N2^2,N2 T ’ * ' “h ^N2 ^^2^8, N 2 Y 

-}- V 02 -^ 02 ^ 1.02 Y 0 2 R 0 2 A 2,02 * * ' T ^ 02 ^ 02 ^ 8502 ^ 

+ ^Ar^Ar Ai 5 Ar ^Ar^Ar-^jAr T * * * + ^Ar^ArAg 5 Ar^ 

+ ^002 ^ 002 ^ 1,002 (^ 2 ) + ^ 002 -^ 002 ^ 2,002 (Jfj + 1 - Y C0 2 R C0 2 A 8,C0 2 t5 ( 5 ) 

By combining like terms in equation ( 5 ), a resultant Cp curve for the gas mixture (air) is generated. 

Cpjair — Ai )a i r (^^2 + ^2, air + ^3, air + + Ay^T 4 + Ag^r^ 5 (6a) 

where 

A^air = -^N2“^N2 ^j,N2 *6 Y Q2 R 02 A jfi2 ^Ar^Ar^,Ar T Y CQ2 R C02 A j,C02 (j 1 to 8 ) ( 66 ) 

In generalized terminology, equations ( 6 a) and ( 6 b) are expressed as 


and 


jmax 

Cp.inix = E A JM* T1 ~ 3 ( 7a ) 

j=l 

n 

^j,mix = ^ 1 YiRiAj i (j = 1, • • • , jmax) (76) 

i— 1 


Gas Expansion Relations 


Parameter 7. With a known curve fit expression for Cp of the gas mixture, the value of 7 for a given 
temperature can be computed directly from its definition and the thermally perfect relation between Cp and 

Cy). 

Cp — Cy — R ( 8 ) 


and 


__ Cp Cp 

7 = — = -~ JL ~ 

°v Cp-R 


( 9 ) 


4 


The gas constant R for the mixture is calculated by dividing the universal gas constant 1 Z by the molecular 
weight of the mixture lT m ; x , 

n 

R = — ( 10 ) 


W„ 


where 


W'rnix - J2X t W t 
1 = 1 

Xi is the mole fraction of each component species determined from 


(ii) 




ft 

A:— 1 * 


( 12 ) 


Parameters u , a, M , and T /Tt- If the flow is assumed to be adiabatic, then the 1-D energy equation 
written for two separate points in the flow field is (ref. 7) 


. u i , u 9 

hi + ~2~ "" + ~ 


(13) 


where the enthalpy h referenced to 0 K is given by 


[ T CpdT 

Jo 


Therefore, 


fTl U ? 2 u, 

Jo CpdT+ i = l 


(14) 


(15) 


If point 2 is selected to represent the stagnation condition, then U 2 — 0 and T 2 = Tt. Then equation (15) 
reduces to 

r T i , 1 ? rTt 

(16) 


f 1 ui r 1 * 

/ Cp dT +-£■=! CpdT 

Jo £ Jo 


or 


rTt 


fTl 


ui f lx 

CpdT - CpdT — CpdT 

J Jo Jo Jti 

where for the selected fifth-order curve fit for Cp (eqs. (7a) and (7b) with jmax = 8) 

fT t 

IT 1 


(17) 


/ 

JT\ 


CpdT = -A X (1 - ^-) + ^2ln ~ + A 3 (T t - Tx)+ ^(t 2 - ??)+ ^(r 3 - T 3 ) 
+ ^ (t 4 - T 4 ) + ^ (T? - Tj 5 ) + ^ (7? - If) 


(18) 


For a specified total temperature, the value of equation (18) can be computed for a range of static temperatures. 
Each value of static temperature T\ (T\ less than Tt) represents a unique point in the expansion of the gas 
from its stagnation conditions. In the temperature increment mode of the TPG program, the user selects the 
total temperature, the desired static temperature increment, and the upper and lower limits for the static 
temperature range. The program simply increments downward through the selected static temperature range 
and computes all of the compressible flow properties for each value of static temperature. The value of T\/Tt is 
therefore not actually computed from any equation but rather is a result of this incremental procedure. From 
equation (17), the expression for the velocity becomes 


u x = 



(19) 
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The speed of sound for a thermally perfect gets is given by 



«i = V'hKTi 

(20) 

and from the definition of Mach number 

Mi — — 

(21) 


a l 


Parameter p/pt» The first law of thermodynamics 



Sq + 6w = de 

(22) 

and the definition of entropy 

Sq Tey — 2 ' ds 

(23) 

along with the thermodynamic relation 

6 

II 

1 

(24) 

yield 

T ds = de + p dv 

(25) 

From the definition of enthalpy, 

h — e 4- pv 

(26) 

differentiation gives 

dh = de + pdv + v dp 

(27) 

The combination of equations (25) and (27) produces 



T ds = dh — v dp 

(28) 

For a thermally perfect gas, 

ii 

^3 

(29) 

The substitution of equation (29) into equation (28) and the definition of specific volume result in 



_ dT dp 

ds = Cn — — 

^ T pT 

(30) 


Substitution of the equation of state p = pRT , which is valid for a thermally perfect gas, gives 


dT dp 

^■s — Cp — - R— (31) 

1 p 

Integration of equation (31) between two different thermodynamic states yields 

*2 — «1 = Y dT ~ K ln W ) (32) 

For an isentropic path from state 1 to state 2, S 2 — s\ = 0 and therefore 



If state 2 represents the stagnation condition, then equation (33) can be rewritten as 


Pi ex p(^ Jti dT ) 


6 


and for the selected fifth-order curve fit for Cp, the closed-form solution to the integral in equation (34) is 

rT, 

>Ti T 


f 


C ?dT = -^ 


^ Jf) ■ M “ 1 \) + ln T + MTt “ Tl)+ If ( T ? ~ T 0 

+ Y C 7 ** - T ?) + T ( T t - T t) + Y ( T ' 5 - r 0 


(35) 


Equations (34) and (35) are used in the TPG code to compute the value of pi/pt for each value of T\. 


Parameter q/pt * A relation for the ratio of dynamic to total pressure q/pt for a given T\ results from the 
following expression 

£ = (36) 

Pt Pl\PtJ 

Then q\/p\ for a thermally perfect gas can be expressed with the use of the equation of state, the definition of 
Mach number, and the equation for the speed of sound as 


Ql _ \p\V\ 
Pi Pi 


i *l 

RT i 




RI\ 


iM, 2 7l 9 


RT 1 

rJui 


(37) 


Parameters p/ pt and A/ A*. With the knowledge of p/pt and T/Tt, the value 
from the equation of state 

of p/pt can be determined 

Pl = piRTi 

(38) 

and 


Pt = PtRTt 

(39) 

Equation (38) is divided by equation (39) to give 


p±_ — ( 

Pt \Pt )\T\ ) 

(40) 


Now the value of A\/A * is obtained by the use of the 1-D continuity equation. 

= p*u*A* 

Equation (41) can be rewritten as 


ulAl =*•■(%:) U - A - 


For isentropic flow, the value of pt is constant (p^i = p^*) and thus equation (42) becomes 

(«)“ iAi= (s) u -' 4 - 


The solution for Ai/A* gives 


A\ _ \pt ) 

^• = te) 


u* 

Ul 


(41) 


(42) 


(43) 


(44) 


As equation (44) indicates, before the value of A\/A * can be computed for any given value of Tj, the values of 
p*/pt and u, must be known. The TPG code is written to calculate all of the compressible flow properties at 
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the sonic point temperature /, before the calculation of these properties at any other T\ values. The value of 
T* is obtained from equations (19), (20), and (21) for the condition M\ - 1 (ui = ai). This gives 

2 f ‘ CpdT = j*RT* (45) 

and equation (9) then yields 

(46) 

Equation (46) can be expressed in terms of only one unknown T* by the use of equations (7a) and (18). An 
iterative technique is implemented in the TPG code to solve for T* in equation (46). 


Normal Shock Relations 

Parameters T 2 , 72 , « 2 > an ^ M 2 . The continuity and momentum equations for 1-D flow across a 
normal shock wave are given in equations (47) and (48), respectively. 

p\u\ = P 2 U 2 (47) 

Pi + Pl^i = P2 + P 2 U 2 ( 48 ) 

The momentum equation is divided by the continuity equation to give 


Pi 

P\u\ 


+ Ul 


P2 

P2 U 2 


+ U2 


(49) 


The relation for the speed of sound, equation (20), can be written with the use of the equation of state as 


a 


2 


TP 

P 


(50) 


The solution of equation (50) for p and substitution into equation (49) yield 

-^L + U1 = -^L +U2 (51) 

Tltil 72«2 

Equation (51) provides the needed relationship across the shock wave. For a given value of 7\, the left-hand 
side of equation (51) is obtained from equations (9), (19), and (20). The right-hand side of equation (51) 
apparently has three unknowns 72 , 112 , and < 22 , but in reality these three variables are all functions of only T 2 . 
(Actually, u 2 depends on both T 2 and T t> 2 , but because a shock wave is considered adiabatic, ht,2 = 1 ; this 

translates to T t p = T ti 1 for a thermally perfect gas because Cp is a function of temperature only.) Therefore, 
from equations (7a), (9), (18), (19), and (20), the right-hand side of equation (51) can be expressed as a single 
elementary, nonlinear function of TV An iterative technique has been implemented in the TPG code to solve 
equation (51) for T 2 . With the known value of T 2 , equations (9), (19), (20), and (21) determine the values of 
72 , U 2 , « 2 , and M 2 , respectively. The desired shock relations, T-j/T) and u 2 /ui , can then be computed. 

Parameters P 2 /P 1 and. P 2 I P\- The 1-D momentum equation in the form 

P2 - Pi = Pl u l - P2«2 ( 52 ) 


and the continuity equation expressed as 


P 2 = 


pim 

U2 


give 


P2 ~Pl = P\u\{u 1 -u 2 )= Piu? 



(53) 

(54) 
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( 55 ) 


Equation (54) is divided by p\ ) and from equation (50) p\/p\ is expressed as 7i/(ai) 2 which gives 

Pi V UlJ 

The expression for p^/pi is then 

Pi V «1 / 

where the values of 71 and M\ are obtained from equations (9) and (21), respectively, and the value for u 2 /ui 
is obtained from equation (51) as discussed in the previous section. With a known value of u 2 /ui, the value of 
P2/P1 is determined directly from the continuity equation. 


P2 _ u\ 
Pi u 2 


(57) 


Parameters Pt^/Pt,l an d Pl/Pt,2 * The value of the total pressure loss Pt^/Pt, 1 can now be computed 
from a combination of static and total pressure ratios. 


Pt, 2 
PU 



(58) 


The values of pt,i/p\ and Pt$/P2 are determined from equation (34) evaluated at T\ and T 2 , respectively, while 
P2/P1 is obtained from equation (56). Similarly, the pitot-static tube pressure ratio p\/pt ,2 is calculated from 


Pi 

Pt2 



(59) 


TPG Code Description 

An interactive computer code, herein referred to 
as the TPG code, has been written based on the 
equations described previously. This code delivers 
a complete table of results within seconds when run 
on a minimally configured workstation. The purpose 
of the code is the output of tables of compressible 
flow properties for a thermally perfect gas or mixture 
of gases styled after the tables of NACA Rep. 1135. 
Output may be in the form of a single table of gas 
expansion properties, a single table of gas expansion 
properties and normal shock properties, or separate 
tables of gas expansion and normal shock properties. 
The entries within the tables are based upon con- 
stant decrements of static temperature or constant 
increments of Mach number. Such entries in terms 
of static temperature decrements reveal a fundamen- 
tal difference between the data of the TPG code and 
NACA Rep. 1135. This fundamental difference is 
highlighted when the TPG code output tables, dis- 
played as constant decrements of static temperature, 
are compared with the calorically perfect gas table 
of NACA Rep. 1135, a single table of compressible 


flow properties for air as a function of Mach num- 
ber. The properties of thermally perfect gases vary 
with both total temperature X) and local static tem- 
perature T rather than with only the ratio of T/7). 
Thus, the utility of the code is its capability to gen- 
erate tables of compressible flow properties for any 
total conditions for any specified range of static tem- 
peratures (T < 7 ) , of course) or Mach numbers. In 
this section, the segments of the TPG code, Ver- 
sion 2.2, are described, code execution procedures 
are included, and examples of all user-supplied in- 
puts are provided. The various types of output are 
also described. 

The FORTRAN code, listed in appendix A, con- 
sists of a Main program, three include files, several 
subroutines, and a Block Data initialization rou- 
tine and is designed to run on a minimally config- 
ured computer workstation. The subroutines arc' 
grouped according to purpose within several files and 
are called as needed from various locations within 
the Main program as well as from other subroutines. 
Table 1 summarizes the files and subroutines which 
comprise the TPG code. A Unix Makefile has been 
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Table 1. TPG Code Version 2.2 Files and Subroutines 


Files 

Contents 

TPG.f 

Version.h 

params.h 

dimens, h 

trangej.f 

cpsubs.f 

ntgrat.f 

titer. f 

Cratio.f 

blockd.f 

initd.f 

Main program 

Include code version identifier 
Include parameters 
Include common blocks 
Subroutine trangej 

Subroutines cpeval, cpNtgra, and cptNtgr 
Subroutines intCp, and intCpT 
Subroutines tliter, and t2iter 
Subroutine Cratio 
Block Data 
Subroutine initd 


developed for compilation of the code and creation 
of an executable file. In the Unix environment, stan- 
dard output prompts the user for the required input 
parameters for code execution, which are written to 
standard input. With one exception, all interactive 
inputs occur in the Main program. Defaults are pro- 
vided for all inputs. A complete execution requires, 
at a minimum, an answer to each interactive prompt 
of a comma followed by a carriage return. Such a 
minimal execution generates a single output file con- 
taining a table of the basic ga s expansion properties 
for the standard composition of air. Generation of 
more complete tables or files for postprocessing re- 
quires specific inputs, which will be described in the 
following sections. 

Main Program 

The first part of the Main program contains the 
interactive inquiries and responses which specify the 
gas components, the desired total temperature, the 
desired outputs, and the output formats. A typical 
example of code execution is shown in table 2 in 
which the user responses are italicized. A response of 
a comma indicates acceptance of a program default. 
Each user response is followed by a carriage return 
(not shown). 

The first response required is an identifier charac- 
ter string (ID), or symbol name, for inclusion within 
the output files. The ID can have a maximum of 
70 characters all entered on a single line; an ID that 
includes blank spaces must be enclosed in quotes. A 
default ID for air may be specified by a comma or 
the word DEFAULT (all uppercase or all lowercase let- 
ters); in this case, the next three questions are omit- 
ted and defaults are assumed. A particular database 
file may be specified next which describes the nec- 
essary thermochemical data for the gases of interest. 
The file name can have a maximum of 40 charac- 


ters and quotes are not required on most comput- 
ers. The contents of such a database file will be 
described in a later section. The code contains the 
database necessary to define a four-species mixture 
of air; this default is chosen by entering a comma. 
The gas mixture definition is completed by the third 
and fourth responses in which the number of individ- 
ual species and the corresponding mass fractions are 
specified. Defaults may be selected at any time, and 
correspond to the standard composition of air. Re- 
gardless of the number of fractions, the code checks 
to see if the sum of the mass fractions equals unity 
and if the database contains sufficient information 
for the requested case. If the mass fractions do not 
sum to unity within a small tolerance, the code is- 
sues an error message and stops execution. If the 
sum is within unity plus or minus the tolerance, the 
last species fraction is recomputed so that the total 
is precisely unity and a warning message is output 
which describes the change. If the database contains 
insufficient information, an error message is output 
and execution is terminated. 

The first two responses are sensitive to the case of 
the alphabetic character entered. That is, upper and 
lowercase alphabetic characters have unique mean- 
ings. However, all succeeding inputs are insensitive 
to case and either upper or lowercase responses have 
the same meaning within the code. 

The next inputs define the composition and for- 
mat of the table(s). The Total Temperature specifies 
the maximum possible upper temperature limit (for 
which the Mach number is zero). Following that, 
is an input that specifies whether the tables are to 
be in terms of static Temperature T or Mach num- 
ber M increments. The default entries in the ta- 
ble(s) are given for incremental Mach numbers from 
an Initial Mach up to an Upper Mach number limit. 
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Table 2. Sample TPG Code Interactive Inputs 


Thermally Perfect Gas Properties Code 
TPG, Version 2.2 


Enter the ID ('character string’) for the gas mixture: 

[default = 4-species Air] 

‘Air: containing standard mass fractions of N2, 02, Argon , and C02’ 

Enter the name of the Database File for the Species properties: 

7 

Enter the number of species composing the mixture: [default = 4] 

7 

Enter the mass fraction, Y, of all 4 component species 
in the order of the Database File (sum=1.0): 

[defaults = 0.7553, 0.2314, 0.0129, 0.0004] 

7 7 7 7 

Enter the Total Temperature (in degrees K) : [default = 400.] 

1500 . 

Increment the table entries by constant Temp (K) or constant Mach number? (T/M) 
[default = M] : 

t 

Enter the Initial Temperature, the Static Temperature Increment, 
and the Low Temperature Limit (in degrees K) : 

[defaults = 1500., 10.0, 10.0] 

1000., 20., 50. 

Do you want to calculate the Normal-Shock Properties? (y/n) : [default = N] 

y 

Enter the number of Tables to be output, 1 or 2 : [default = 1] 


Enter the Tabular output location: screen=6, file=9 [default = 9 : TPG. out] 

7 

Do you want to output a TECPL0T file [TPGpost.dat] ? (y/n): [default = N] 

y 

Do you want to Ratio the data to the case of a 
Calorically Perfect Gas ? (y/n) : [default = N] 

y 

Enter the Calorically-Perf ect Constant Gamma: [def ault=l . 4] 

1.4 
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Table 3. Sample TPG.out Header Information 


| Thermally Perfect Gas Properties Code | 
j TPG, Version 2.2 | 

Table of Thermally Perfect Gas Flow Properties 

Gaseous Mixture: Air: 4-Species Mixture of N2, 02, Argon, and C02 

Database file name: air. data 

Species Names are: 

N2 , 02 , Ar , C02 , 

Species Mass Fractions are: 

0.75530 0.23140 0.01290 0.00040 
Species Mole Fractions are: 

0.78092 0.20946 0.00935 0.00026 
Mixture Properties: Molecular Wt . - 28.9663 

Gas Constant = 2.87035E+02 J/(kg*K) 

Stagnation Temperature = 1500.000 (K) 


Alternatively, if the constant Temperature option is 
specified, the table entries are for Static Tempera- 
ture Increments from an Initial Temperature down 
to the Low Temperature Limit. The initial values 
are the total conditions by default. These inputs are 
separated by commas and any, or all, may be de- 
faulted simply by entering a comma. Gas expansion 
properties are always included in the table and nor- 
mal shock properties may be obtained by an answer y 
for “yes” to the inquiry. The relative locations of the 
expansion and shock property tables are determined 
by the answer to the next inquiry with regard to the 
number of tables (either 1 or 2). The normal shock 
properties may be output in the same table as the 
expansion properties in a wide format (158 charac- 
ter width), or may be output separately, in a second 
table cross-referenced to the first table of expansion 
properties by the Mach number and static temper- 
ature. The table of normal shock properties begins 
with sonic conditions and includes all subsequent (su- 
personic) entries. Because the tables can be rather 
extensive, they are written by default to an output 
file TPG.out. However, in response to the ques- 
tion of tabular output location, the tables could be 
written to standard output, which is typically the 
computer screen. 

The next inquiries determine the output of addi- 
tional files for data postprocessing and plotting. Be- 
cause of extensive usage at NASA Langley Research 
Center, the Tecplot 1 (ref. 8) code format has been 


^Tecplot: trademark of Amtec Engineering, Inc., 

Bellevue, WA 98009-3633. 


chosen for these files. The format of this ASCII file, 
which is described in reference 8 and in appendix B, 
could be easily modified to satisfy the input require- 
ments of many different graphics postprocessors. A 
single zone of information for all variables and all ta- 
ble entries is included. For subsonic conditions the 
normal shock properties have no physical meaning 
and are defaulted to a large negative number as 
a warning flag. The first postprocessor file TPG- 
post.dat contains the same basic data as that in the 
TPG.out tabular output. The second postprocessor 
file TPGratio.dat is generated only if a comparison 
of the thermally perfect properties to calorically per- 
fect properties is desired. In this file all data are 
output as a ratio of thermally perfect values normal- 
ized by the corresponding calorically perfect values. 
A specific constant 7, which is required for the calor- 
ically perfect calculations, is the response to the fi- 
nal prompt. This last interactive input is the only 
one not requested from the Main or at the beginning 
of the interactive execution. The input is requested 
when and if the normalization calculation is about to 
begin. 

Following the interactive inputs, the database 
file, if one is requested, is opened, read, and closed. 
Certain properties of the gas mixture are computed 
and the coefficients of the c p relationship for the 
mixture are defined (eq. (7a)). The TPG.out file 
is opened and header information is written to it 
stating the gas mixture specifications and the total 
temperature. Sample header information is shown 
in table 3. Initialization of the case is completed 
by calculation of the sonic conditions as required by 
equation (44). 
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A program loop for the temperature or Mach 
range between the desired limits computes the ba- 
sic gas expansion properties. These properties are 
then output unless the single table format with the 
normal shock properties has been requested. A sec- 
ond loop for all supersonic Mach numbers calculates 
these normal shock properties. After the second loop 
is completed, the tabular output to TP G. out is com- 
pleted and the file is closed. The same information is 
then output in postprocessor format to the TPG- 
post.dat file. Finally, normalization of the ther- 
mally perfect data with that of a calorically perfect 
gas is performed and the normalized results are out- 
put to the TPGratio.dat file. Execution of the code 
is completed by closing all files and recording warning 
summaries of the times the valid temperature ranges 
have been exceeded during the calculations. These 
warnings, if any, are related to the database infor- 
mation and are written to an additional file called 
Extrap* WARN for later reference. 

Include (or Header) Files 

The Version, h file contains a character variable 
which specifies the current version number of the 
computer code. This information is output at the 
beginning of each interactive execution and to the 
first lines of the TPG.out header. 

Two files params.h and dimens. h contain in- 
formation which determines the size of the computer 
memory required to execute the TPG code and, thus, 
the size of the case which may be specified. The di- 
mensions of the primary arrays within the code as 
defined below are included in a FORTRAN param- 
eter statement in params.h and have the following 
values as released in Version 2.2: 

1. mtrm (=8) number of polynomial 

terms which define 
relationship for Cp/ R 

2. mspc (=10) maximum number of gas 

species allowed 

3. mevf (=3) maximum number of 

temperature ranges 
per species allowed and 
corresponding sets of 
polynomial coefficients 
of equation (2) 

4. mtner (=10000) maximum number 

of computed static 
temperatures 

All values assigned to the preceding parameters, 
except mtmi, may be changed at anytime to suit 


particular user requirements. The code must then 
simply be recompiled and linked to form a new ex- 
ecutable file. A change in mtrm may require code 
changes in the subroutines found in cpsubs.f which 
perform specific calculations using the c v (T)/R re- 
lationships. Changes to mtrm which are permitted 
without recoding the routines in cpsubs.f are de- 
scribed in “Polynomial Summation and Integration 
Subroutines. ” 

The dimens. h file specifies common blocks con- 
taining all of the major arrays that are dependent on 
the parameters described above. Additional common 
blocks include the data for sonic conditions, certain 
species and mixture information, and several charac- 
ter variables. 

FORTRAN include statements allow these 
header files to be separated from each FORTRAN 
routine so that the definitions are included at com- 
pile time as if the individual statements were within 
each routine. For example, the code dimensions need 
to be specified only once rather than in each routine. 

Temperature Rangefinder Subroutine 

The specification of any single function which can 
accurately describe the relationship between heat ca- 
pacity and static temperature over a wide tempera- 
ture range for typical gases is an extremely difficult 
task. However, when limited ranges are considered, 
maintenance of continuity and smoothness of sepa- 
rate functions at the limits of the individual ranges 
are much simpler. References 5 and 6 provide polyno- 
mial curve fits of the specific heat, specific enthalpy 
per degree, and specific entropy for a number of gases 
for multiple temperature ranges. The older report 
(ref. 6) includes a large selection of gases defined by 
five-term, fourth-order polynomials and listed as ac- 
curate for temperature ranges of 300 to 1000 K and 
1000 to 5000 K. Continuity of the polynomials for 
each thermodynamic quantity was maintained at the 
common point of 1000 K. Appendix 2 of reference 5 
contains updated polynomial expressions for a subset 
of the gases of reference 6. These updated polynomi- 
als contain eight terms and include all integer powers 
of T from —2 to 5. The published temperature ranges 
for these polynomials are 200 to 1000 K and 1000 to 
6000 K. Detailed information regarding the deriva- 
tion and validity of these functions may be found in 
references 5 and 6. 

To evaluate equation (7a) for specific heat at 
a given temperature, the appropriate temperature 
range is required; the trangej subroutine was de- 
veloped to provide this. Given a particular temper- 
ature, the routine searches to determine the correct 
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temperature range from the set of valid ranges for 
the gas mixture. For a temperature beyond the lim- 
its of any range, the closest range is specified for 
extrapolation and a warning message is output to 
the file Extrap. WARN for that particular tempera- 
ture. A single extrapolation warning message is writ- 
ten at the end of the TPG.out file which notes the 
temperature below which extrapolations were per- 
formed. This warning message appears when any 
of the component species of a gas mixture requires 
extrapolation. 

Polynomial Summation and Integration 

Subroutines 

The equations defined in the derivation section 
require evaluation of Cp, the integral of Cp, and the 
integral of c p /T. These evaluations are performed 
numerically in three subroutines contained in file 
cpsubs.f. These three subroutines contain code 
specific to the form of the polynomial expression of 
equation (2): the cpeval subroutine evaluates Cp 

as a function of T from equation (7a), the cpNtgra 
subroutine evaluates the integral of Cp for the range 
of T to Tt from equation (18), and the cptNtgr 
subroutine evaluates the integral of Cp/T for the 
range of T to Tt from equation (35). Required 
inputs to the integration subroutines are the mixture 
polynomial coefficients, the static temperature, and 
the total temperature. 

For integrations that cover multiple temperature 
ranges, cpNtgra and cptNtgr must be called once 
for each temperature range and the corresponding 
set of coefficients. This function is performed in 
the TPG code by two subroutines intCp and intCpT 
in the file ntgrat.f. Individual results of calls to 
the cpNtgra and cptNtgr subroutines are summed 
to yield a total integral. The trangej subroutine 
determines the ranges which include T and Tt. The 
total integral is divided into separate integrals as 
follows: 

fTt fT hm rTt 

/ CpdT= / CpdT + / CpdT (60) 
JT Jt JT hm 

This decomposition of the single integral is then 
performed as required to cover all of the temperature 
ranges. 

The cpeval, cpNtgra, and cptNtgr subroutines 
assume a polynomial expression of the form given 
in equation (7a). The lowest power of T is -2. 
The highest power of T must be at least 2, which 
requires the parameter mtrm > 5. Any mtrm > 5 
may be set in params.h with no changes required 
in the FORTRAN code. However, for mtrm < 5 or a 


different expression for c p /R , these three subroutines 
must be modified. An increase of mtrm beyond the 
code release value of 8 would include higher order 
powers of T in the polynomials. 

Nonlinear Iteration Subroutines 

As stated previously, the properties at sonic con- 
ditions must be known in order to compute the value 
A\/A* for any T\. A combination of equations (19), 
(20), and (21) gives the relationship for T\ as 

2 / cp dT = A/jSi RTi (6 1 ) 

M 

for any M\ , in which 71 is a function of T\ . At sonic 
conditions M\ = M* = 1 by definition and equa- 
tion (61) becomes equation (45). An application of 
Newton’s method (ref. 9) to equation (61) yields an 
iterative technique for T\ that has been implemented 
in subroutine t liter (file titer. f). The derivative in 
Newton’s method is defined by a first-order finite dif- 
ference expression rather than by an analytic expres- 
sion. For the iterative solution of this well-behaved 
function, an initial guess of Ti = 0.8*T t was satisfac- 
tory. The algorithm is called once in the TPG code to 
determine the static temperature at sonic conditions; 
however, the algorithm is actually valid for determin- 
ing the static temperature for any Mach number, not 
just M*. Therefore, t liter is also called for each 
Mach number tabular entry when the Mach number 
increment option is selected. 

Subroutine t2iter also implements a Newton it- 
eration to solve equation (51) for the static temper- 
ature T2 behind a normal shock. This function is 
not as well-behaved as that for T\ and tends to be- 
come rather flat near the minimum. The solution 
by Newton’s method either converges slowly or be- 
comes unstable; for some conditions, desirable levels 
of convergence cannot be attained. Checks on the 
magnitude of the derivative and the number of itera- 
tions have been implemented in the code to prevent 
divergence; the desired convergence tolerance is in- 
creased, if needed, to prevent divergence. A warning 
message is written to standard output each time this 
tolerance is changed. For maximum stability, the it- 
eration is started with an initial guess defined by the 
minimum of either T2 = 0.99 * Tt or the calorically 
perfect value based upon local conditions. 

One of the arguments for each Newton subroutine 
is the desired convergence tolerance epsil. The 
tolerance is set in a Main program data statement 
at an acceptable value for a computer word size of 
32 bits. For computers that use a different word size 
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(e.g., a CRAY 2 with a 64-bit word), this tolerance 
may be changed for improved accuracy or stability. 

Calorically Perfect Comparison 

Subroutine Cratio is called after all gas proper- 
ties have been computed from the thermally perfect 
relationships. With the last user input 7 , this subrou- 
tine applies the appropriate equations from NACA 
Rep. 1135 to compute corresponding gas properties 
for a calorically perfect gas. Each property com- 
puted by the thermally perfect equations is normal- 
ized by the calorically perfect value to allow analysis 
of the magnitude of the thermal dependency. Ra- 
tios of physically meaningless values, such as nor- 
mal shock properties at subsonic values of M\ y are 
set to defaults of large negative numbers in Ver- 
sion 2.2 of the code. The ratios are output to the 
file TPGratio.dat for postprocessing in the Tecplot 
format described in appendix B. Again, this format 
can be easily changed to suit the needs of a differ- 
ent postprocessor. All file formats are included in 
Cratio. 

Data Initialization 

Data initialization is located in three places 
within the TPG code. A Block Data subroutine con- 
tains data statements which initialize the variables 
associated with individual species and the gas mix- 
ture. Most of the Block Data is associated with the 
default database information for the standard com- 
position of air. A second set of data statements are 
included within the Main program itself. These state- 
ments set defaults for the variables that do not need 
to be passed to subroutines via common blocks such 
as the universal gas constant, the sonic Mach num- 
ber, the Newton iteration convergence tolerance, and 
a finite representation of infinity for numeric output. 
Additionally, the default responses for several of the 
interactive inputs (particularly the default yes or no 
answers) are set in the Main data statements. 

The third location for data initialization is in the 
subroutine initd where the arrays of gas properties 
are set to defaults. The TPG code could be modified 
to allow the user to redefine the total temperature, 
the temperature or Mach increments, and/or the 
temperature or Mach limits, and then recompute the 
tables for a given mixture within a single execution. 
For such a modification, all the necessary arrays 
could be reset to their defaults by a single subroutine 
call. 
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Database File Format 

The thermochemical data required by the TPG 
code for a given mixture of gases are defined in a 
database file to be read at execution time. As pre- 
viously stated, a database for the standard compo- 
sition of air is contained within the code and may 
be accessed as the default. However, for mixtures of 
other gases or more complete models of air, a sep- 
arate database file may be provided. The format 
is simple; a two- line data file header is followed by 
groups of data for each species. A sample two-species 
database is shown in table 4, and specific database 
files are included with the discussions of the test cases 
to follow. 

The first line is a descriptor of the data to follow 
on the second line and, as such, is not actually read 
by the TPG code. Those data consist of the number 
of species for which the file includes thermochemi- 
cal data and the number of temperature ranges for 
which separate polynomails define c p / R as a function 
of T. These integer values are read in free format and 
are separated by one or more blanks. Version 2.2 of 
the code requires that, for each species within a par- 
ticular file, each of the polynomials must have the 
same limiting temperatures for each range, except for 
the absolute minimum and maximum temperatures 
of the overall definition. That is, if the polynomials 
for one gas are valid from 200 to 1000 K, and from 
1000 to 5000 K, then each gas within that database 
file must also be defined by one polynomial up to 
1000 K, and a second one above 1000 K. The abso- 
lute extrema of 200 and 5000 K are not required to be 
identical since the code simply extrapolates the poly- 
nomials beyond these limits. Extrapolation warnings 
are based upon the worst case extrema of the most 
conservatively defined gas. The requirement that all 
of the temperature ranges coincide is not a limita- 
tion of the basic algorithm but only of Version 2.2 
of the code. (Note that this limitation can be eas- 
ily overcome by splitting any particular temperature 
range into multiple consecutive ranges with identical 
coefficients.) 

The next set of lines are repeated for each of the 
gaseous species in the file; the number of individual 
gaseous species expected by the code is given in the 
second line. Error handling has been incorporated 
within the code to recognize End-of-File on most 
computers without terminating execution, (i.e., the 
case in which the actual number of species defined 
is less than expected). However, some computers 
and/or FORTRAN compilers may not function prop- 
erly if the header information and the actual number 
of gases do not agree. 
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Table 4. Sample Thermochemical Database File 


# of species # of Temperature ranges: NASP TM 1107 + Mods. 

2 2 
Name: Molecular Wt . 

N2 28.0160 

tmin tmax 
20. 1000. 

Cp/R Coef f icients : cl(-2) — >cl(5) 

- 1 . 33984200E+0 1 1 . 34280300E+00 3 . 45742000E+00 5 . 74727600E-04 

-3 . 2171 1900E-06 7 . 50775400E-09 -5 . 90150500E-12 1 . 50979900E-15 
tmin tmax 
1000. 6000. 

Cp/R Coefficients : c2(-2) — >c2(5) 

5 . 87702841E+05 -2 . 23921563E+03 6 . 0668697 1E+00 -6 . 13957913E-04 
1 . 49178026E-07 -1 . 92307130E-11 1 . 06193594E-15 0 . 00000000E+00 

Name: Molecular Wt . 

02 32.0000 

tmin tmax 
30. 1000. 

Cp/R Coefficients : cl (-2) — >cl(5) 

3 . 88517500E+01 -2 . 70630800E+00 3 . 56119600E+00 -3.32782400E-04 

-1 . 18148000E-06 1 . 10853500E-08 -1 .49299400E-11 5. 99553800E-15 

tmin tmax 
1000. 6000. 

Cp/R Coef f icients :c2(-2) — >c2(5) 

- 1 . 05642070E+06 2 . 41 123849E+03 1 . 73474238E+00 1 . 31512292E-03 

-2.29995151E-07 2 . 13144378E-11 -7 . 87498771E-16 0 . 00000000E+00 


The thermochemical data for the first gas begins 
with the third line of the file, which is the descrip- 
tor for the next line of data. The fourth line con- 
tains the name of the gas (five characters), one or 
more blanks, and its molecular weight. Each poly- 
nomial definition follows beginning with the lowest 
temperature range. Line six specifies the minimum 
and maximum temperatures for the first polynomial. 
After another descriptor line, the coefficients A\ to 
As (or Amtrm if mtrm ^ 8) for the polynomial are 
tabulated, in free format separated by blanks or on 
multiple lines. After a descriptor line, the minimum 
and maximum temperatures for the second polyno- 
mial are specified next, followed by a descriptor line 
and a second set of tabulated coefficients. This pat- 
tern is repeated for all the temperature ranges spec- 
ified in the database header. The data for a second 
gas follow directly after these lines in the same for- 
mat as the first and continue in this pattern until all 
gases have been defined. 

The TPG code will use the gas definitions in 
the order specified within the database file and the 
user must input mass fractions in that same order. 
However, mass fractions of zero are acceptable inputs 


if one or more of the leading species within a given file 
are to be omitted. Also, the polynomial coefficients 
for a particular gas do not all have to be nonzero. In 
particular, a calorically perfect gas may be defined 
by a database in which all of the coefficients are 
zero except A3, which equals the caloric constant for 

Cp/R- 

TPG Code Validation 

Calorically Perfect Temperature Regime 

The first validation test of the TPG code was the 
verification of accuracy in the calorically perfect tem- 
perature regime in which air is the test gas. In this 
temperature regime, the specific heat of air is nearly 
constant and therefore the TPG code results should 
be nearly identical to the tables of NACA Rep. 1135. 
For this test case, the default data file for standard 
four-species air was used (table 5) along with the 
default (standard) values for the mass fraction com- 
position of air. The code default total temperature 
of 400 K, which is well within the calorically per- 
fect temperature regime for air, was selected and the 
case was run with the TPG code in the tempera- 
ture increment mode. The upper temperature limit 
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(Initial Temperature) was defaulted to the total tem- 
perature. The Static Temperature Increment and 
Low Temperature Limit were chosen to be 2 K and 
20 K, respectively. (The Static Temperature Incre- 
ment merely determined the number of rows of tab- 
ulated output in the TPG.out file. The choice for 
a Low Temperature Limit determined the maximum 
Mach number presented in this file.) The TPG.out 
file for this test case is given in table 6 (a). Table 6 (b) 
shows the results for the same test case obtained with 
the Mach number increment option of the TPG code. 
A Mach number increment of 0.05, a minimum Mach 
number of 0.0 (the default value), and a maximum 
Mach number of 10.0 were selected for this case. A 
comparison of table 6 (b) with tables I and II from 
NACA Rep. 1135 for any matching Mach number 
shows the desired agreement and validates the TPG 
code in this calorically perfect temperature regime. 
The small differences that are observed between the 
corresponding values of table 6 (b) and the tables of 
NACA Rep. 1135 are in either the third or fourth 
digit. To better examine these small differences, the 
TPG code generated the ratio of the gas properties 
calculated with the thermally perfect relations (ta- 
bles 6 (a) and 6 (b)) to those calculated from the calor- 
ically perfect relations (tables of NACA Rep. 1135). 
These results have been plotted versus Mach number 
in figures 1 and 2. From these figures, the maximum 
difference between the values listed in table 6 (b) and 
the tables of NACA Rep. 1 135 appears to be approxi- 
mately 0.5 percent. Note that these small differences 
actually represent the small amount of error associ- 
ated with the tabulated values of NACA Rep. 1135 
which resulted from the calorically perfect gas as- 
sumption. Close inspection of table 6 (a) or table 6 (b) 
reveals that a slight variation with temperature ac- 
tually exists in the specific heat ratio, which leads to 
the noted differences when compared with the calor- 
ically perfect gas properties. 

Thermally Perfect Temperature Regime 

The first step in verification of the TPG code 
in the thermally perfect temperature regime was to 
compare the 7 variation for air computed by the 
TPG code with both the 7 variation calculated from 
NACA Rep. 1135 (eq. (180)) and the theoretical ther- 
mally perfect 7 variation from reference 4 obtained 
with the theoretical c p data for N 2 , O 2 . Ar, and CO 2 . 
The variation of 7 as a function of temperature is 
presented in figure 3. Below 300 K, 7 values predicted 
by NACA Rep. 1135 have attained their asymptotic 
value of 1.400 while the 7 values predicted by the 
TPG code have reached a slightly greater asymptotic 
value of 1.401, which is in agreement with the theo- 
retical data. Above 300 K, both the TPG code and 


the exponential expression for 7 in NACA Rep. 1135 
produce nearly identical 7 variation that is in excel- 
lent agreement with the theoretical 7 values up to 
approximately 1600 K; at this point, the 7 curve of 
NACA Rep. 1135 begins a gradual divergence away 
from the TPG code curve of 7 and the theoretical 
results. This divergence results from the inability of 
the exponential expression of NACA Rep. 1135 (es- 
sentially a single constant 0 exponential curve fit) to 
accurately model the true thermally perfect behavior 
of the gas at the higher temperatures. The TPG code 
avoids this limitation by utilization of multiple-range 
polynomial curve fits to model the variation of spe- 
cific heats; this particular application of the code con- 
sisted of two eight-term, fifth-order polynomial curve 
fits. From a true gas behavior standpoint, neither 
NACA Rep. 1135 nor the TPG code predicts an ac- 
curate 7 variation for air above 1600 K because of the 
effects of dissociation. A discussion of this tempera- 
ture limitation associated with the application of the 
TPG code is presented in the next section. 

The remaining gas properties computed by the 
TPG code in the thermally perfect temperature 
regime were verified for test cases at total tempera- 
tures of 556 K (1000° R), 1111 K (2000° R), 1667 K 
(3000° R), and 2778 K (5000° R) for air (standard 
four-species composition) as the test gas. These four 
specific test cases were selected so that a comparison 
could be made with the caloric imperfection charts 
of NACA Rep. 1 135 which give the compressible flow 
gas properties for air computed from the thermally 
perfect relations based upon equation ( 1 ) of this pa- 
per. Each chart from NACA Rep. 1135 shows the 
variation of one gas property, as normalized by its 
calorically perfect value, versus Mach number for the 
above mentioned four total temperatures. The TPG 
code results for these four test cases were also nor- 
malized by their calorically perfect values and then 
plotted together with the results from the caloric 
imperfection charts of NACA Rep. 1135. (The re- 
sults of NACA Rep. 1135 have been regenerated here 
for plotting purposes from the appropriate equations 
in ref. 1.) Figures 4 15 present the results of the 
TPG code and NACA Rep. 1135 and correspond to 
charts 9 20, respectively, from reference 1. Exami- 
nation of figures 4 15 shows that, overall agreement 
(maximum difference less than one percent) exists for 
the compressible flow properties computed from the 
two different methods. The small differences that are 
observed in these figures can be attributed in most 
cases to one of the two sources of 7 inaccuracies noted 
previously. The predominate difference noted in fig- 
ures 4 15 is between the curves of the TPG code 
and NACA Rep. 1135 for the test case at a total 


17 



temperature of 2778 K (5000° R). This difference is 
discernible from the start of the expansion (clearly 
evident in figs. 4, 8, and 12) and continues through- 
out the Mach number range. This results from the 
difference in the predicted 7 values between the two 
methods for the temperature above approximately 
1600 K. (See fig. 3.) The other noticeable differ- 
ence between the curves of the TPG code and NACA 
Rep. 1135 occurs at the total temperatures of 1667 K, 
1111 K, and 556 K. At some given value of Mach 
number, the curves of the TPG code diverge from 
the asymptotic plateau value of NACA Rep. 1135. 
This divergence is most apparent in figures 6-9, and 
15. The onset of divergence begins at approximately 
M = 3 for T t = 556 K, at M = 5 for T t = 1111 K, 
and at M = 6 for Tt = 1667 K. Examination of the 
TPG. out files for these test cases (not included in 
this report) revealed that the Mach numbers at this 
onset of divergence for the above three test cases all 
correspond to a static temperature value of approx- 
imately 200 K. The divergence of the curves of the 
TPG code and NACA Rep. 1135, which is discernible 
in the plotted data at about 200 K, is actually a con- 
sequence of the difference in the predicted 7 values as 
shown in the low-temperature range of figure 3. This 
comparison verifies the TPG code accuracy with the 
derived thermally perfect gas relations based upon 
polynomial expressions for c p . Although these test 
cases were all for air, which is a primarily diatomic 
gas, the TPG code is also valid for polyatomic gases. 
The utilization of a polynomial curve fit for c p makes 
the TPG code applicable to any molecular structure 
of the gas. 

As a verification of the above statement, figure 16 
presents the 7 variation for a triatoinic molecule 
CO2 predicted by the TPG code versus the theoret- 
ically thermally perfect 7 variation from reference 4. 
Unlike figure 3, there is no 7 curve from NACA 
Rep. 1135 because its equation (180) is not valid for 
a triatomic molecule such as CO2. The agreement in 
the 7 variation in figure 16 is excellent, which is just 
a confirmation of the curve fit for c p as a function of 
the temperature utilized by the TPG code. Table 7 
gives the TPG. out file for this CO2 test case. (Note 
that for the CO2 test case, the air. data database file 
was used with the mass fractions for N2, O2, and Ar 
set equal to zero.) 

Temperature Boundaries for Application 

of TPG Code 

The TPG code will provide valid results as long 
as its application is kept within the thermally per- 
fect temperature regime for the gas of interest. If 
the code is used outside of the region for which c p 


is a function of temperature only, then the ther- 
mally perfect results produced by this code will no 
longer accurately reflect what actually occurs in na- 
ture. To illustrate what happens outside the bound- 
aries of the thermally perfect temperature regime, 
comparisons have been made with air data from 
NBS “Tables of Thermal Properties of Gases.” (See 
ref 3.) Figure 17 presents the variation of specific 
heat (nondimensionalized by the gas constant) with 
temperature for air up to 3000 K. This figure con- 
tains data from the TPG code, NACA Rep. 1135 
(eq. (178)), and reference 3. The results of the TPG 
code and NACA Rep. 1135 are nearly identical; the 
NBS data (plotted for three pressures of 0.1, 1.0, and 
10 atm) coincide with the other results only between 
approximately 500 K and 1500 K. Above 1500 K, the 
NBS data yields significantly larger values of Cp/R 
and clearly becomes a function of pressure as well as 
temperature. This phenomenon is a result of dissoci- 
ation effects. Below 500 K, the NBS curve for 10 atm 
diverges to greater values of Cp/R . In this region of 
relatively low temperatures and high pressures, real 
gas effects are the cause of this deviation from the 
thermally perfect behavior. (Here the term real gas 
refers to a gas for which the intermolecular forces 
are not negligible.) From reference 3, examination of 
the compressibility factor table for air in this range 
of temperature and pressure verifies this conclusion. 
(Note that at low temperatures and also relatively 
low pressures, the thermally perfect results will ac- 
curately predict the actual behavior of the gas prop- 
erties.) The heat capacity information given in fig- 
ure 17 is presented in an alternate form in figure 18, 
which shows the variation of the ratio of specific heats 
with temperature for air from the same three sources. 
The characteristics discussed above are also observ- 
able in this figure. These two figures emphasize that 
the TPG code user must remain aware of the low- 
and high-temperature boundaries associated with the 
thermally perfect assumption for the particular gas 
under consideration. (Note that the user should not 
confuse the above issue with the extrapolation warn- 
ings which result from exceeding the limits of appli- 
cation of a particular curve fit for c p /R as a function 
ofT.) 

Additional Sample Cases 

Although only air and CO2 test cases have been 
presented thus far, the TPG code has the potential 
to produce tables of compressible flow properties for 
any gas or any mixture of gases for which the required 
polynomial curve fit information for Cp is available. 
Four additional sample cases have been included in 
this section to provide a better understanding of the 
full capability of the TPG code. The four gases are 
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steam, a mixture of steam and air, a mixture of 
tetrafluoromethane and argon, and nitrogen, which 
are examples of monatomic, diatomic, triatomic, and 
polyatomic gases. The database files for these four 
sample cases are given in tables 8-11. Tables 12-15 
are the TPG.out files which contain the tabulated 
output of compressible flow gas properties for these 
four sample cases. The values of mass fractions and 
the stagnation temperatures were arbitrarily selected 
for all four cases. Stagnation temperature values 
were checked to ensure that they were safely below 
the temperature limits for all component species 
at which dissociation effects would invalidate the 
thermally perfect assumption. The low-temperature 
limit in all cases was simply selected as the tmin 
value of the lowest temperature range curve fit for 
c v /R. For the multiple gas species sample test 
cases, the highest low-temperature range tmin values 
were selected as the low-temperature limit. This 
ensured that no curve fit extrapolations' would be 
encountered. 

A verification that the results from the low- 
temperature region of tables 12-15 are not in the 
region of real gas effects cannot be performed until a 
value of total pressure is specified. This value deter- 
mines the magnitude of the static pressures from the 
TPG.out file. In the case of a mixture of gases, the 
value of static pressure for the mixture needs to be di- 
vided into the values of partial pressures for the com- 
ponent gas species by application of the mole fraction 
information printed in the header of the TPG.out 
file. From the combination of static temperature and 
pressure for each component gas species, the com- 
pressibility factor Z is obtained from a source such 
as reference 3, which provides tables of compressibil- 
ity factors for several gas species. If the value for 
Z differs by more than one percent from unity, then 
real gas effects should be considered if accurate re- 
sults are desired. A check to determine if liquefaction 
of the gas has occurred for the given static temper- 
ature and pressure is also appropriate at this stage 
in the analysis. The possibility of liquefaction ex- 
ists when the compressibility factor table does not 
provide data at the desired ranges of low tempera- 
ture and high pressure. A table or plot of the vapor 
pressure versus temperature will verify whether the 
particular component species is in either the liquid 
or gas phase. 


Conclusions 

A set of compressible flow relations for a ther- 
mally perfect gas has been derived for a value of Cp ex- 
pressed as a polynomial function of temperature and 
developed into a computer program referred to as the 
TPG code. The code produces tables of compress- 
ible flow properties similar to those found in NACA 
Rep. 1135. Unlike the tables of NACA Rep. 1135 that 
are valid only in the calorically perfect temperature 
regime, the TPG code results are also valid in the 
thermally perfect temperature regime, which gives 
the TPG code a considerably wider range of temper- 
ature application. Accuracy of the TPG code in the 
calorically perfect temperature regime was verified 
by comparisons with the tables of NACA Rep. 1135. 
The small differences observed in this comparison 
were identified as an inaccuracy of the tabulated 
values in NACA Rep. 1135 that resulted from the 
calorically perfect gas assumption. In the thermally 
perfect temperature regime, the TPG code was ver- 
ified by comparisons with results obtained from the 
method of NACA Rep. 1135 for calculation of ther- 
mally perfect compressible flow properties. The TPG 
code is applicable to any type of gas and is not re- 
stricted to only diatomic gases as is the method of 
NACA Rep. 1135. In addition, the TPG code can 
accept any specified gas mixture for which the com- 
ponent gas species polynomial curve fit information 
for c p is known because the calculation of the perti- 
nent thermochemical mixture properties is performed 
within the code. 

The temperature limits for application of the 
TPG code have also been examined. The thermally 
perfect behavior of any gas species will terminate at 
some combination of low temperature and high pres- 
sure due to the presence of real gas effects and at 
some value of high temperature due to dissociation 
(except for monatomic gases). Beyond these tem- 
perature limits, the TPG code may no longer accu- 
rately predict the actual behavior of gas properties. 
However, within the temperature range in which the 
thermally perfect assumption is valid, the TPG code 
provides accurate gas species or gas mixture com- 
pressible flow properties. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
May 3, 1994 
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c write(12 , 9300) ( (A { 1 , j ) , 1 = 1 , mt nn) , j = l, ncvf ) avel (k) = sqrt (gamma (k) *r*t (k) ) 

C9300 format (lp2elS. 5) m(k) * zero 

c cpti(k) = zero 

c prat (k) = 1 . /exp(cpti (k) /r) 

c Opening an output file & Writing Input Summary rhorat(k) = prat (k) /trat (k) 
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subroutine intCp ( t , tstag, tmin, tmax, A, ncvf , cpin) subroutine cpeval {A,t,cp) 

C23456789 123456709 123456789 123456789 123456709 123456789 123456789 12 
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Subroutine for the calculation of the integral of cp/T subroutine trangej (t, tmin,tmax, jmax, j) 

from T to Tstag : C23456709 123456789 123456789 123456709 123456789 123456789 123456709 12 
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Alternate Subroutine Entry point for output of Summary 
Information describing Total Extrapolation Warnings. 
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Write Long Output: 

wr itelll, 9 02 4) mst , Rgamm, Rprat , Rrrat , Rtrat , Rqopt , Raast , Rvast , 
Rm2 , Rp2pl , Rr2rl , Rt2tl , Rpt 21 . Rpltl 
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mzts*0 .0/ 

data tmin, tmax/0 .0, mcvf*0.0/ 



Appendix B 

Tecplot ASCII Output File Formats 

The TPG code output data files for postprocess- 
ing are in the Tecplot ASCII format defined in ref- 
erence 8 and summarized below. This format was 
chosen for its simplicity and current widespread us- 
age at NASA Langley Research Center. The Tecplot 
interactive plotting program, developed by Amtec 
Engineering, Inc., provides a suitable means for vi- 
sualization of the tabular results of the TPG code. 
However, the output file formats could be easily mod- 
ified for compatibility with many other graphics post- 
processing software packages; and the ASCII output 
format facilitates the transfer of data between var- 
ious computer platforms. The following discussion 
will be useful in the conversion process. 

The TPG output files which utilize this format 
are the TPGpost.dat and the TPGratio.dat files. 
TPGpost.dat contains the basic properties of the 
thermally perfect gas, which are also output to the 
TPG.out file. A sample TPGpost.dat file is shown 
in table Bl. The TPGratio.dat file contains the 
same basic properties but in the form of ratios with 
the properties of a corresponding calorically perfect 
gas. The ratio of specific heats of the calorically 
perfect gas is one of the interactive inputs at the time 
of TPG code execution. 

The first line in each of these postprocessor data 
files is a TITLE record containing the ID for the cor- 
responding TPG calculation. A record designation 
TITLE = precedes the data title; a sample TITLE 
record for the default data is shown in the table. The 
ID is entered as part of the interactive input and is 
the portion within the double quotes. 

The second line in the postprocessor files contains 
the names of the variables to be included, and also 
the number of variables for which data are included 
in the file. A record designator of VARIABLES = 
precedes the list of names. The current version of 
Tecplot restricts the length of each name to five char- 
acters. The names in the TPGpost.dat file (and 
the data to follow) are listed in the same order as 
found in the basic TPG.out tabular listing. The 
TPGratio.dat file utilizes the same name s and or- 
der, except that the variables T\ and (3 = — 1| 

are omitted. Because the caloric data are not a func- 
tion of temperature, the Mach number is identical 
for both calorically and thermally perfect flows. The 
variable names are abbreviated forms of the head- 
ing names for the columns of tabulated data listed 


in TPG.out. For example, all variables which are 
ratios have the / omitted from the name. Pressure, 
temperature, and density (RHO) variables utilize sin- 
gle uppercase letters to designate the function (i.e., 
P, T, and R). The dynamic pressure is represented 
by q. A lowercase t represents the subscript t for 
total, or stagnation, conditions. The conditions up- 
stream and downstream of a shock are signified by 
the numerals 1 and 2, respectively; the 1 is occa- 
sionally omitted for brevity. The letters str for * 
denote sonic conditions. The same variable names 
are used in both TPGpost.dat and TPGratio.dat; 
however, the values within TPGratio.dat are ratios 
with respect to calorically perfect gas calculations. 
The VARIABLE record in the TPGpost.dat file thus 
appears as shown in the second line in the table. If 
the normal shock properties are not computed, the 
last six variables are omitted. 

The third line of the files (table Bl) is referred 
to as a ZONE record. A ten character zone title, 
identified by T=, follows the record designator ZONE. 
The TPG zone title contains the total temperature 
for which the data were computed. The number of 
data points 1= recorded within the file follow the 
zone title. Although the Tecplot format allows IJ- 
ordered data sets, the TPG code output requires only 
the single I index. Therefore, the next data output 
is fixed for all cases as J=l. The last information on 
the ZONE record describes the order in which the data 
that follow are listed. The F=P0INT signifies that the 
data are presented in variable sequence first. That 
is, for each Mach number and/or temperature data 
point, the N variables are listed sequentially, followed 
by the N variables of the next Mach number and/or 
temperature data point. 

The actual gas property data follow the three 
header lines. The data are written in FORTRAN 
exponential format with six significant digits; the 
page format provides up to a maximum of eight 
values per line. Thus, when normal shock data 
is included, the sixteen values require exactly two 
successive lines in TPGpost.dat. Successive Mach 
number and/or temperature data points always begin 
on new lines. Data values are output for every 
variable at each data point. Data values for the 
normal shock relations at subsonic Mach numbers 
(which have no physical meaning) are defaulted to 
a large negative number as a warning flag. This 
flag is intended to help the user recognize and avoid 
plotting these results. The flag magnitudes are set 
in the initd subroutine for TPGpost.dat and in 
subroutine Cratio for TPGratio.dat. 


32 


Table Bl. Sample TPGpost.dat File 


TITLE = "Air: 4-Species Mixture of N2, 02, Argon, and C02" 

VARIABLES = Ml , T1 , Gamma, PIPt, RIRt , TITt , Beta, qlPt, AAstr, VAstr, M2 , P2P1, R2R1 , T2T1, Pt2Pt, 


PlPt2 

ZONE T= "Tt = 2000.". 1= 6 


0 . OOOOOE+OO 
9 . 99990E+07 
5 . 93104E-01 
1 . 20218E+00 
8 . 59855E-01 
1 . 01864E+00 
1. OOOOOE+OO 
1. OOOOOE+OO 
1 . 08104E+00 
1 . 00554E+00 
1 . 28333E+00 
1 . 06318E+00 


2 . 00000E+03 
0. OOOOOE+OO 
1 . 90000E+03 
6. 19233E-01 
1 . 80000E+03 
8.74541E-01 
1 . 73804E+03 
1 . OOOOOE+OO 
1 . 70000E+03 
1 . 06954E+00 
1 . 60000E+03 
1 . 23309E+00 


J= 1, F=POINT 

1 . 29801E+00 1. OOOOOE+OO 1. OOOOOE+OO 1. OOOOOE+OO 1. OOOOOE+OO 0 . OOOOOE+OO 
-9 . 99990E+07 -9.99990E+07 -9.99990E+07 -9.99990E+07 -9.99990E+07 -9.99990E+07 
1 . 30003E+00 8 . 00246E-01 8.42364E-01 9.50000E-01 8.05126E-01 1.82981E-01 
-9 . 99990E+07 -9.99990E+07 -9.99990E+07 -9.99990E+07 -9.99990E+07 -9.99990E+07 
1 . 30225E+00 6 . 33527E-01 7.03919E-01 9.00000E-01 5.10538E-01 3.04987E-01 
-9 . 99990E+07 -9.99990E+07 -9.99990E+07 -9.99990E+07 -9.99990E+07 -9.99990E+07 
1 . 30376E+00 5.44935E-01 6.27069E-01 8.69019E-01 0. OOOOOE+OO 3.55233E-01 
1. OOOOOE+OO 1 . OOOOOE+OO 1. OOOOOE+OO 1. OOOOOE+OO 1 . OOOOOE+OO 5.44935E-01 
1 . 30474E+00 4 . 95617E-01 5.83079E-01 8.50000E-01 4.10660E-01 3.77853E-01 
9 . 26497E-01 1.19181E+00 1 . 14390E+00 1.04189E+00 9.99393E-01 4.95919E-01 
1 . 30754E+00 3.82658E-01 4.78322E-01 8.00000E-01 8.04318E-01 4.12012E-01 
7.91595E-01 1 . 73683E+00 1.52015E+00 1 . 14254E+00 9.81353E-01 3.89929E-01 
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Table 5. Database File for Standard Four-Species Composition of Air 


# of species #of Temperature ranges : NASP TM 

4 2 

Name: Molecular Wt . 

N2 28.0160 

tmin tmax 

20. 1000. 

Cp/R Coefficients: cl(-2) --> cl(5) 

-1. 33984200 E+ 01 1 . 342803 00E+00 3 . 45742000e+00 5, 

-3 . 217 11900E-06 7 . 50775400E-09 - 5 . 9 0150500E- 12 1, 

tmin tmax 

1000. 6000. 

Cp/R Coefficients: c2(-2) --> c2(5) 

5 . 877 0284 1E+ 05 -2 . 23921563E+03 6 . 06686971E+00 -6, 

1 . 49 17802 6E-07 - 1 . 923 07 13 0E- 11 1 . 06193594E-15 0. 

Name: Molecular Wt . 

02 32.0000 

tmin tmax 

30. 1000. 

Cp/R Coefficients: cl(-2) --> cl(5) 

3 . 88517 500E+ 0 1 -2 . 70630800E+00 3 . 56119600E+00 -3, 

-1.18148000E-06 1 . 10853 500E- 08 -1 . 49299400E-1 1 5. 

tmin tmax 

1000. 6000. 

Cp/R Coefficients: c2{-2) --> c2{5) 

-1. 0564207 0E+ 06 2 . 4 1 123 849E+03 1 . 73474238E+00 1. 

-2 .2999 515 IE -07 2 . 13 144378E-11 -7 . 8749877 IE-16 0. 

Name: Molecular Wt . 

Ar 39.9440 

tmin tmax 

20. 1000. 

Cp/R Coefficients: cl(-2) --> cl(5) 

O.OOOOOOOOE+OO O.OOOOOOOOE+OO 2 . 50000000E+00 0. 

O.OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO 0. 
tmin tmax 

1000. 6000. 

Cp/R Coefficients: c2(-2) --> c2(5) 

O.OOOOOOOOE+OO O.OOOOOOOOE+OO 2 . 50000000E+00 0. 

O.OOOOOOOOE+OO O.OOOOOOOOE+OO 0.00000000E+00 0. 

Name: Molecular Wt . 

C02 44.0220 

tmin tmax 

80. 1000. 

Cp/R Coefficients: cl(-2) --> cl(5) 

-5. 34531900 E+ 03 2 . 287 854 0 0E+02 2 . 56174500E-01 1. 

-2 . 09949800E-05 1 . 40430800E-08 -3 . 81732100E-12 0. 

tmin tmax 

1000. 6000. 

Cp/R Coefficients: c2(-2) --> c2(5) 

1 . 154 60 08 1E+0 5 -1 . 78337370E+03 8 . 28644239E+00 -8. 

4 . 26107946E-09 -1 . 81443266E-12 6 . 29130739E-16 0. 


1107 + Mods . 


74727600E-04 
509799 0 0E- 15 


13 9579 13E-04 
OOOOOOOOE+OO 


3278240 OE-O 4 
99553 800E- 1 5 


31512292E-03 

OOOOOOOOE+OO 


OOOOOOOOE+OO 

00000000E+00 


OOOOOOOOE+OO 

OOOOOOOOE+OO 


68213800E-02 

00000000E+00 


983 5694 5E-05 
OOOOOOOOE+OO 



Table 6. TPG.out File for Calorically Perfect Temperature Regime Test Case 
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320.0 1.3991 4.562E-01 5.702E-01 0.8000 0.5072 4.012E 



Table 6(a). Continued 
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(b) TPG code run in Mach number increment mode 
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Table 8. Database File for Steam Sample Case 


# of species lof Temperature ranges 

1 2 
Name: Molecular Wt . 

H20 18.01528 

tmin tmax 

80. 1000. 

Cp/R Coefficients: cl(-2) --> c 1 { 5 ) 

8 . 08886500E+02 -3 . 43014000E+01 4 . 51171200E+00 

8. 64754500E-06 -7 . 33842100E-09 2 . 39461300E- 12 

tmin tmax 

1000. 6000. 

Cp/R Coefficients: c2(-2) --> c 2 { 5 ) 

1 .0349719SE+06 - 2 . 41269814E+03 4 . 64611030E+00 

-6.8368312SE-07 9 . 42 64 7001 E- 1 1 -4 . 82238112E-15 


-3 . 21581800E-03 
O.OOOOOOOOE+OO 


2 . 29199858E-03 
O.OOOOOOOOE+OO 


Table 9. Database File for Mixture of Steam and Air Sample Case 


# of species *of Temperature ranges 

5 2 

Name: Molecular Wt . 

H20 18.01529 

tmin tmax 

80. 1000. 

Cp/R Coefficients: cl(-2) --> cl (5) 

8 . 08886500E+02 -3 . 43014000E+01 4 . 51171200E+00 

8 . 64754500E-06 -7 . 33842100E-09 2 . 39461300E-12 

tmin tmax 

1000. 6000. 

Cp/R Coefficients: c2 ( - 2 ) --> c2(5) 

1 . 03497195E+06 -2 . 41269814E+03 4 . 64611030E+00 

-6. 83683125E-07 9 . 42 647001E- 1 1 -4 . 82238112E-15 

Name; Molecular Wt . 

N2 28.0160 

tmin tmax 

20. 1000. 

Cp/R Coefficients: cl ( -2 ) — > cl(5) 

-1 . 33984200E+01 1 . 34280300E+00 3 . 45742000e+00 

-3 . 21711900E-06 7 . 50775400E-09 -5 . 90150500E-12 

tmin tmax 

1000. 6000. 

Cp/R Coefficients: c 2 ( - 2 ) — > c2(5) 

5 . 87702841E+05 -2 . 23921563E+03 6 . 06686971E+00 

1.49178026E-07 -1 . 92307130E-11 1 . 06193594E-15 

Name: Molecular Wt . 

02 32.0000 

tmin tmax 

30. 1000. 

Cp/R Coefficients: cl(-2) --> cl (5) 

3.88517500E+01 -2 . 70630800E+00 3 . S6119600E+00 

-1.18148000E-06 1 . 10853500E-08 -1 . 49299400E-11 

tmin tmax 

1000. 6000. 

Cp/R Coefficients: c2(-2) --> c2(5) 

-1 .05642070E+06 2 . 41123849E+03 1 . 73474238E+00 

-2 . 29995151E-07 2 . 13144378E-11 -7 . 87490771E-16 

Name: Molecular Wt . 

Ar 39.9440 

tmin tmax 

20. 1000. 

Cp/R Coefficients: cl(-2) — > cl (5) 

O.OOOOOOOOE+OO O.OOOOOOOOE+OO 2 . 50000000E+00 
O.OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO 
tmin tmax 

1000. 6000. 

Cp/R Coefficients: c2(-2) — > c2(5) 

O.OOOOOOOOE+OO O.OOOOOOOOE+OO 2 . 50000000E+00 
O.OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO 
Name: Molecular Wt . 

C02 44.0220 

tmin tmax 

80. 1000. 

Cp/R Coefficients: cl ( -2 ) — > cl (5} 

-5 . 34531900E+03 2 . 28785400E+02 2 . 56174500E-01 

-2 . 09949800E-05 1 . 40430800E-08 -3 . 81732100E-12 

tmin tmax 

1000. 6000. 

Cp/R Coefficients; c2(-2) — > c2(5) 

1 . 15460081E+05 -1 . 78337370E+03 8 . 28644239E+00 

4 . 26107946E-09 -1 . 01443266E-12 6 . 29130739E-16 


-3 . 21 581800E-03 
0.00000000E+00 


2 . 29199858E-03 
O.OOOOOOOOE+OO 


5 . 74727600E-04 
1.50979900E-15 


-6.139S7913E-04 

O.OOOOOOOOE+OO 


-3 . 32782400E-04 
5.99553800E-15 


1 . 31512292E-03 
O.OOOOOOOOE+OO 


O.OOOOOOOOE+OO 

O.OOOOOOOOE+OO 


0.00000000E+00 

0.00000000E+00 


1 . 68213800E-02 
O.OOOOOOOOE+OO 


-8 . 98 3 5 6945 E- 05 
O.OOOOOOOOE+OO 


Table 10. Database File for Mixture of Tetrafluoromethane and Argon Sample Case 

# of species #of Temperature ranges 

2 2 
Name: Molecular Wt . 

CF4 88.0110 

tmin tmax 

75. 1000. 

Cp/R Coefficients: cl(-2> — > cl (5), Mod 5th Order 

0 . 000000e+00 2.187449E+01 2.959023E+00 5.807661E-03 

6 . 277310E-05 -1 . 517045E-07 1.350004E-10 >4 . 297511E- 14 

tmin tmax 

1000. 6000. 

Cp/R Coefficients: c2 ( -2 ) --> c2(5). Mod 5th Order ( JANAF) 

0 . 000000e+00 -2 . 441536e+03 1.500538e+01 -8 . 5433 55e-04 

1 . 942 454e-07 -2 . 232766e-ll 1.016662e-15 0.000000e+00 

Name: Molecular Wt . 

Ar 39.9440 

tmin tmax 

20 . 1000 . 

Cp/R Coefficients: cl{-2) -- > cl(5) 

O.OOOOOOOOE+OO O.OOOOOOOOE+OO 2 . 50000000E+00 O.OOOOOOOOE+OO 
0 .OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO 
tmin tmax 

1000. 6000. 

Cp/R Coefficients: c2(-2) --> c2(5) 

O.OOOOOOOOE+OO O.OOOOOOOOE+OO 2 . 50000000E+00 0.00000000E+00 

O.OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO O.OOOOOOOOE+OO 


Table 11. Database File for Nitrogen Sample Case 

* of species #of Temperature ranges 

1 2 
Name: Molecular Wt . 

N2 28.0160 

tmin tmax 

20 . 1000 . 

Cp/R Coefficients: c 1 ( - 2 ) --> c 1 ( 5 ) 

-1. 33984200 E+ 01 1 . 34 280300E+00 3 . 45742000e+00 5 . 74727600E-04 

-3 . 2171 1900E-06 7 . 50775400E-09 -5 . 90150500E-12 1 . 50979900E-15 

tmin tmax 

1000. 6000. 

Cp/R Coefficients: c2(-2) --> c 2 ( 5 > 

5 . 077O2841E+O5 -2 . 23921563E+03 6 . 0668 6971E+00 - 6 . 13957913 E- 04 

1 .49178026E-07 -1 . 92 307130E-11 1 . 06193594E-15 0.00000000E+00 
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Table 15. TPG.out File for Nitrogen Sample Case 
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Figure 4. Effect of caloric imperfections on ratio of local speed to speed of sound at M = 1. 
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Figure 9. Effect of caloric imperfections on ratio of local cross-sectional area of stream tube to cross-sectional 
area at M — 1. 
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Figure 10. Effect of caloric imperfections on static temperature ratio across normal shock wave. 
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Figure 11. Effect of caloric imperfections on static density ratio across normal shock wave. 
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Figure 12. Effect of caloric imperfections on rati 
pressure downstream. 
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Figure 13. Effect of caloric imperfections on static pressure ratio across normal shock wave. 
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